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Abstract: We consider the leptonic and semileptonic (SL) lepton flavour violating (LFV) 
,-£h ' decays of the charged leptons in the minimal supersymmetric standard model (MSSM). 

The formalism for evaluation of branching fractions for the SL LFV charged-lepton decays 
with one or two pseudoscalar mesons, or one vector meson in the final state, is given. 
Previous amplitudes for the SL LFV charged-lepton decays in MSSM are improved, for 
instance the 7-penguin amplitude is corrected to assure the gauge invariance. The decays 
are studied not only in the model-independent formulation of the theory in the frame of 
MSSM, but also within the frame of the minimal supersymmetric SO (10) model within 
which the parameters of the MSSM are determined. The latter model gives predictions 
for the neutrino-Dirac Yukawa coupling matrix, once free parameters in the model are 
appropriately fixed to accommodate the recent neutrino oscillation data. Using this un- 
ambiguous neutrino-Dirac Yukawa couplings, we calculate the LFV leptonic and SL decay 
processes assuming the minimal supergravity scenario. A very detailed numerical analysis 
is done to constrain the MSSM parameters. Numerical results for SL LFV processes are 
given, for instance for r~ — ► e~ (//" ) ir° , r~ — ► e~(fi~)r}, t~ — ► e~([i~)r}', t~ — ► e~(fi~) p°, 
t~ — ► e~ (/x~ ) <fi, t~ — > e~ ) lj, etc. 
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1. Introduction 

The neutrino oscillation experiments gave a first experimental evidence beyond the stan- 
dard model (SM) of the electroweak interactions. In the SM, neutrinos are massless, purely 
left-handed particles, so there is no leptonic analogy of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix in the SM. The neutrino oscillation experiments proved that the neutrinos 
do mix and that they do have mass. The mixing matrix in the lepton sector, the Maki- 
Nakagawa-Sakata (MNS) matrix has bi-large mixing structure, indicating that the source 
of the lepton-flavour mixing is different from the corresponding mixing in the quark sector. 
The lepton-flavour mixing observed in neutrino oscillation is the first confirmation that 
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the lepton flavour is not conserved quantity. Therefore, experimental observation of the 
other lepton- flavour violating (LFV) processes is naturally expected. Theoretical study 
of such processes has a long history before the observation of neutrino oscillations. The 
model independent study of the operators using the SM fields [g, g, § shows that there 
are no LFV operators of dimension less or equal to four. There is one dimension five 
LFV operator that induces neutrino oscillations. The LFV decays can be induced only 
with the operators of dimension six or more. As the new physics is expected to appear at 
the scale much larger than the electroweak scale ~ 246 GeV, the LFV decay effects are 
expected to be much more suppressed than the neutrino oscillation effects. Model inde- 
pendent study of the LFV processes gives the limits on LFV which every model has to 
satisfy. A model dependent analysis depends on the structure of the model, but are much 
more predictive than the corresponding model independent analysis. Therefore, the both 
approaches are indispensable for a theoretical study of LFV. Although the leptonic LFV 
processes have been studied extensively both in model independent way and using various 
models [J§ §, @, the semileptinic (SL) LFV processes have been studied only in few models 

The super symmetric (SUSY) models have much nicer theoretical properties than their 
non-SUSY counterparts. For example, quadratically divergent contributions to the Higgs 
boson mass from heavy (e.g. GUT scale) particles cancel with its SUSY partners and as 
a result the gauge-hierarchy problem is much better resolved. The supersymmetrization 
of the SM cannot be done without additional assumptions. For instance, in the super- 
symmetric version of the SM there are dimension four operators violating both the lepton 
number (^) and the baryon number leading to very fast proton decay. That led us to 
the introduction of a discrete Z2 symmetry, so called, i?-parity to forbid such undesirable 
terms. The SUSY breaking also has to be done in such a way not to induce too large 
flavour violation effects. There are few successful SUSY breaking mediation mechanisms, 
such as gravity mediation [14], gauge mediation ]i~5| , anomaly mediation [jOI], gaugino me- 
diation [fT^], radion mediation ||18U , etc. The best established among them is the minimal 
supergravity model (mSUGRA) [14|] which assumes that the SUSY breaking occurs in the 
hidden sector at very high scale, which communicates with the visible sector (containing 
SM) with flavour-blind gravitational interactions. The induced soft SUSY breaking mass 
terms are requested to be universal at the SUSY breaking mediation scale (say, the (re- 
duced) Planck scale), and are therefore flavour-diagonal. The magnitude of the soft SUSY 
breaking mass terms obtained is in the range to induce potentially observable consequences 
in the visible sector. The renormalization group (RG) flow from the (reduced) Planck scale 
to the mass scale of the right-handed neutrinos, induces the flavour non-diagonal terms 
in the SUSY soft breaking terms for the sleptons, through the flavor- non diagonal Dirac- 
neutrino Yukawa matrices they contain |19[ . They can lead to considerable LFV effects, 
which, depending on the model parameters can be in the range of the forthcoming LFV 
experiments [2C]. 

In this paper, we assume the MSSM with three right-handed neutrinos as the low- 
energy effective theory below the GUT scale. In such a framework, neutrino oscillation 
data suggest the existence of very massive right-handed neutrinos which give rise to the 
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small left-handed neutrino masses through the see-saw mechanism [21]. In SO(W) models, 
the required right-handed neutrinos can be naturally embedded into the common multiplet 
together with the SM particles for each generation. In this paper, the minimal renormazible 
SUSY 50(10) model |24|, |25|, |26|, 27] will be taken as a starting theoretical frame. 

One of the advantageous points of this model is the automatically conserved i?-parity 
defined as R = (— i~j3{B-L)+2S ^ w \ iere g represents the spin of a field. Namely, the 50(10) 
model discussed here spontaneously breaks the gauged B — L symmetry by two units, 
leading to an automatic i?-parity conservation. The breaking of SO(10) group to the SM 
gauge group, SU(3)c x SU(2)l x U(1)y [^g, 29] and its phenomenological consequences 
H 31] has already been discussed in our previous publications. 

The main goal of this paper is an analysis of neutrinoless SL LFV decays of charged lep- 
tons within the MSSM model where parameters are obtained from the underlying SO (10) 
model. At the same time we intend to see how the previous phenomenological analyses 
constrain the LFV parameters. The paper consists basically of three parts which are given 
in three sections. In section two, we give the MSSM form factors comprised in the LFV 
amplitudes at the quark-lepton level. We rederive these form factors, because some of them 
were not derived completely in the previous literature. In section three, the charged-lepton 
SL LFV amplitudes at lepton-meson level are derived using a simple hadronization proce- 
dure for the quark currents. The branching ratios corresponding to these amplitudes are 
also given. In section four, the minimal renormalizable SUSY 5*0(10) model is described, 
the parameters of the MSSM model are derived from the 5O(10) model, too. Using the 
predicted Yukawa couplings in the minimal 5*0(10) model, we perform a numerical estima- 
tion for the SL LFV processes. The last section is devoted to a summary. In Appendix A, 
we give our notation for the mass matrices of the neutralino, chargino and sfermions. The 
MSSM Lagrangian for fermion-sfermion-(gaugino, Higgsino) interaction and the trilinear 
interactions with Z boson are given in Appendix B and C, respectively. In appendix D we 
present the loop functions needed to evaluate the SL LFV processes. The quark content 
of meson states, essential for the hadronization of quark currents, is listed in Appendix 
E, together with constants that define the hadronized quark current in 7-penguin and 
Z-boson-penguin amplitude. 



2. Effective Lagrangian for LFV Interactions 
2.1 Sources for LFV Interactions 

Even though the soft SUSY breaking parameters are flavour blind at the scale of the SUSY 
breaking mediation, the LFV interactions in the model can induce the LFV sources at 
low-energy through the renormalization effects. In the following analyses we assume the 
mSUGRA scenario ]Jl4|] as the SUSY breaking mediation mechanism. At the original scale 
of the SUSY breaking mediation we impose the boundary conditions on the soft SUSY 
breaking parameters, which are characterized by only five parameters: too, M1/2, Aq, B 
and \x. Here, mo is the universal scalar mass, M 1 / 2 is the universal gaugino mass, and Aq is 
the universal coefficient of the trilinear couplings. The parameters in the Higgs potential, B 
and /i, are determined at the electroweak scale so that the Higgs doublets obtain the correct 
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electroweak symmetry breaking VEV's through the radiative breaking scenario [32]. The 
soft SUSY breaking parameters at low energies are obtained through the RGE evolutions 
with the boundary conditions at the GUT scale. 

Although the SUSY breaking mediation scale is normally taken to be the (reduced) 
Planck scale or the string scale (~ 10 18 GeV), in the following calculations we impose the 
boundary conditions at the GUT scale (~ 10 16 GeV), and analyze the RGE evolutions 
from the GUT scale to the electroweak scale. This ansatz is the same as the one in the 
so-called constrained MSSM (CMSSM). 

The effective theory which we analyze below the GUT scale is the MSSM with the 
right-handed neutrinos. The superpotential in the leptonic sector is given by 

W Y = Y^(u R ) i q j H u + Y^(d R ) i q j H d 

+ Y^v^ijHu + YiHeWjHi + \m Rij {u c r Uu r ) + fiH d H u , (2.1) 

where the indices i, j run over three generations, H u and H d denote the up- type and down- 
type MSSM Higgs doublets, respectively, and Mr., is the heavy right-handed Majorana 
neutrino mass matrix. We work in the basis where the charged-lepton Yukawa matrix Y e 
and the mass matrix M Rij are real-positive and diagonal matrices: Yg J = Y ei Sij and M Ri . = 
diag(M Rl , M R2 , M Ra ). Thus, the LFV is originated from the off-diagonal components of 
the neutrino-Dirac- Yukawa coupling matrix Y v . The soft SUSY breaking terms in the 
leptonic sector is described as 

t f^2\ ~. , -t f m 2\ ~ . Tf / 2 



- £ soft = 9i (™|L qj + u Rl (ml) u Rj + d Ri (m|) _ _ d R 



+ A \ m i).. h + D m ( m l)ij "no + 4i ("»!)« £Rj 
+ m 2 Hu HlH u + m 2 Hd H ] d H d + ^B/iH d H u + X -B v M Ri] v R ^ R j + h.c\ 
+ {^AqjHu + A^qjHa + h.c? 



+ (^M\BB + ^M 2 W a W a + ^M 3 G a G a + h.cA . (2.2) 

As discussed above, we impose the universal boundary conditions at the GUT scale such 
that 



(rnf) = {ml).. = (mj).. = m\hj , 



2 2 2 

m H u = m Hd = m , 

A i i = A Q Y$, A i i = A Q Y i j 

A^ = A Y^, A i i=A Y? 



Ah = M 2 = M 3 = Mi /a , (2-3) 



- 4 - 



and evolve the soft SUSY breaking parameters to the electroweak scale according to their 
RGE's [33]. The fi parameter and the B parameter are determined at the electroweak scale 
so as to minimize the Higgs potential, 



N " tan 2 0-1 2 z > 

B» = ~\ (m Hd + m 2 Hu + 2\fi\ 2 ) sin 2/3 ( = sin 2^ . (2.4) 



The LFV sources in the soft SUSY breaking parameters such as the off-diagonal components 
of \jn^\ and A l J are induced by the LFV interactions through the neutrino Dirac Yukawa 

couplings. For example, the LFV effect most directly emerges in the left-handed slepton 
mass matrix through the RGE's such as 



d ( i\ d 
Mr \ m a).. = I m ) 



dfl \ £ Jij dfJL 



MSSM 



+ ^2 (mfyjY, + Y$Y v m\ + 1Y$m\Y v + 2m 2 H Y^Y v + 2A\A v ) (2.5) 

where the first term on the right-hand side denotes the MSSM term with no LFV. In the 
leading-logarithmic approximation, the off-diagonal components (i ^ j) of the left-handed 
slepton mass matrix are estimated as 

where the distinct thresholds for the right-handed Major ana neutrinos are taken into ac- 
count by the matrix = log (^jf^J We can see that the neutrino-Dirac- Yukawa 
coupling matrix plays the crucial role in calculations of the LFV processes. 

2.2 Effective Lagrangian in terms of quark fields and LFV form factors 

In any model containing standard model as the low-energy effective theory, an effective 
Lagrangian for the SL LFV decays of a lepton contains only three terms: the photon- 
penguin, the Z-boson-penguin and the box term, 

iC cS (£i ^£j + q + q') = i£2s + <S + ^eff X - (2-7) 
These terms have the following generic structure, 

i£j s (x) = -ie 2 J d 4 y£j(x) [{-d 2 llx + @ x d Xfi )D(x - y){V^P L + V^Pr) 

+ <Tp,er x D{x - y)Q%p L + v^p r )] ii{x) 

x Y, QM-fQiv), ( 2 - 8 ) 

q=u,d,s 

m z c w 
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q{h q ~ 2Q g s 2 w h» - h q ris)q(x), (2.9) 



x 

q=u,d,s 



q a qb=uu,dd,ss,ds,sd 



+ B^ aqb (£^P L £,)(q altl P R q b ) + B% aqb (£^P R £ t )(q a ^P L q b ) 
+ B^ljPLliKqaPLto) +B* aqb {l j P R £ i ){q a P R q b ) 

+ B^ii^PLli^qaa^PRqb) + B^&a ^P^iWa^ 'P L q b ) 



(2.10) 



where sw = sin 6^ and c\y = cos 6w, Qq is the quark charge in units of e, and I^ q 
is weak quark isospin. D{x — y) is the Green function for the massless scalar particle, 
contained in the photon propagator. The structure of the photon-penguin term in the 
effective Lagrangian is a consequence of the gauge invariance. Especially, the first term 



must contain <f) x d x ^, which was neglected in reference [2C]. The information on the model 
under consideration is contained in the form factors Va-} R , a = 1,2, V^ ,R , B^ aq , and 
B a q aqb , a = 1,2, 3, 4. In the following three subsections these form factors are given for the 
MSSM. 

2.2.1 The photon-penguin form factors 

The amplitude for £{ —* Ijj* for an off-mass-shell photon process is obtained the from 
corresponding part of the effective Lagrangian neglecting the quark current and the photon 
propagator, 



= iT; = -eu e . 



(q 2 ^-q^)(vtPL+v^p R ) + ia^ q »(v£ 'P L +r£ \p r ) 



u tv (2.11) 



The amplitude is written without photon polarization vector. 

In the MSSM the photon-penguin amplitude has two contributions, a chargino and a 
neutralino contribution. That reflects in the structure of the form factors, 

^^W + pW, a = 1,2 (2.12) 

with C and TV" subscript denoting the chargino and neutralino part of a form factor. 

Because of the gauge invariance, the zeroth-order and the first-order term in Tay- 
lor expansion in momenta and masses of incoming and outgoing particles are equal zero. 
Here, the second order term in Taylor expansion is presented, and higher order terms are 
neglected. 

The neutralino contributions are 

1 



~ 576yr 2i i AxI iAX mj 

ex 

o^AxT ' (2 - 13) 

< )B = V™ L \l~r, (2-14) 
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v. 



(N)L 
27 



32vr 2 



<Jjvgw*(-i: 



77? 



ex 



2{x Ax f + 3(*° x ) 2 - 6(x AX ) + 1 - 6(x° x f ln(4 x ; 



+ ly jAX iy iAX ( 



-l)m,- 



6(1 - x° x )4 

1 



2 

m- 

ex 



2(x Ax f + 3(x AX ) 2 - Q(x AX ) + 1 - Q{x Ax f Hx AX ) 



6(1 - x AX r 

1 -( x ax? + ^ + ^ x ax)H x ax) 



m 



(l-x A 



AX) 



V. 



(N)R 
27 



27 



\L*->R, 



(2.15) 
(2.16) 



where x AX = M~ /rrig . The chargino contributions are given by 



V 



(C)L 
1 7 



57 6vr 2 jAX iAX m 



2 



16 - 45(x^ x ) + m{x Ax f - 7(x AX ) 3 + 6(2 - 3(x AX )) ln(x AX ) 



{l-x 



AX J 



V 



v. 



V 



1-y |L<->it) 



(C)-R 
17 

(C)£ _J_ 
" 32vr 2 

2 + 3(x 



^ii(e) r i?(e)* 1 



m 



i'x 



AX> 



6(1 - x AX ) 



, (~<L{e) r L{e)* 



777 



vx 



b AXJ 



6 ( x Ax) + ( x Ax) + Q ( x Ax) H x Ax) 



, pL{e) n R(e)* 
+ U jAX U iAX m y 



6(i-x AX r 

1 -3 + 4(^ x )-(x^) 2 -21n(x^) 



"x 



{1-x 



AX) 



V. 



(C)R 
27 



n 



(C)L, 

2-y |L<-»it) 



(2.17) 
(2.18) 



(2.19) 
(2.20) 



where x 



AX 



M 2 _/mf, . The form factor contributions are written in the same way as 
X A ' Ux y 
in Ref. [ 20 ] , including the loop functions into the expressions, to make the comparison 

with the results of Ref. [20| easy. Both chargino and neutralino parts of the form factors 



agree with the corresponding form factors in the reference [2C] if the terms proportional 
to the mass of the lighter mesons rrij are neglected. Nevertheless, these terms cannot be 
neglected, because the constants N^J^^ and C^ AX ' (see Appendix B) also depend on the 
lepton masses is such a way that in some cases the term proportional to rrij is larger than 
the term proportional to the mass of decaying lepton rrij . 
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2.2.2 The Z-penguin form factors 

The amplitude for l{ — > ijZ* , for the off-mass-shell Z-boson process, (obtained analogously 
as the photon-penguin amplitude) is 



Ml 



(Air) 



,(N)L 



(C)R 



(2.21) 



where j>^ L ' R anc j r p^) L > R are ^ e criar gi no and neutralino part of the total form factors, 



PZ>. Here are the expressions for these form factors: 



V 



{C)L 



V 



n R{e) n R(e)* 
^jBX^iAX 



E BA : "\ r >'\ x 



-Fxfmi 



2 2 * 

m , m / 

Xa Xb ' 



" m \- A ' m h ) + SABG zefi (m: 



2 2 > 



+ \pjBX '- B.\ K i AS 



,R(e) L(x-)r,L(e)* 



(N)L 



R(e) AT R(e) 



-2F 2 (m 



2 2 2 

« y ,m— ,m— 
x A x s 



)]}, 



/V H ; AT 
^jBY ly iAX 



2D 



2 2 < 
ex ) m e> , 



+ S YX G Ze f 2 (m|o ,m. 



2 * 
ex' 



+ 



L(e) ^(x ) N R(e)* 



BX E BA 



AX 



+ N jBX h BA yV --«^ 



V 
V 



(C)R 

z 

(N)R 



vf )L {L 



"iAX 

R), 
R). 



i m x° B m x A Fl< < m 
-2F 2 (m 



2 

ex 1 



m 2 ~n , m-n )1 

X A Xb' 1 



2 2 2 

ix , m -q , m 



fl(x°) 

BA 



-E 



L(X°) 
BA ■ 



(2.22) 



(2.23) 

(2.24) 
(2.25) 



G Ze and G^ e are constants appearing in the SM Ze^ vertices, 



-gi^{G Ze P L + G R e P R } 



'J 



,2 1 

w 



1 1 s 

"9— + 



Pl + 



„2 1 
S W 

CW 



(2.26) 



while ^ and E^'^* ^ are constants in the Z-boson — chargino and Z-boson — neutralino 



(e) 

vertices, and D YX is a constant in the Z-boson — selectron vertex. These constants are de- 
fined in Appendix C. -Fi(a, b, c) and F 2 (a, b, c) are loop functions contained in the triangle- 
diagram part of the amplitude, and f\ and f 2 are the loop functions coming from the 
self-energy part of the amplitude. They are given in Appendix D. 

The terms in Eqs. ( |2.22 ) and Q2.23 ) which have a corresponding contribution in the 
photon amplitude (the leading order photon-penguin amplitude comes from six Feynman 
diagrams, while the Z-boson-penguin amplitude has eight Feynman-diagram contributions) 
have been compared by replacing Z-boson vertices with corresponding photon vertices, and 
an agreement was found. The remaining two Feynman diagram contributions which are 
embraced by curly brackets in Eqs. ( [2.22| ) and ( 2.23 ) have been checked carefully. The 



BA 



new terms in V Z L in comparison with Ref. [^0| are third (self-energy-type term) and 
fourth term. Further, neither of our terms in J>^' R does agree with amplitude in Ref. pp|| , 
although the expression in the curly brackets is almost equal to it. 
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2.2.3 The box form factors 

The box contribution to the SL LFV t — ► (-iq a % amplitude comes from two box-diagrams 
in the leading order of perturbation theory. The box-amplitude reads 



B lq a q b ( U tj I^PlUI, ) (Uq a 7m P £%, ) + B^ (u tj Y ''P R U£. ) ( U q a ^P R V qb ) 



q a qb=uu,dd,ss,ds,sd 

+ B^ a ^(ut j 'fP L Ui i )(ug 9 'y tl P I(0qi ) + B^ aqb {u e ^P R u ti )(u qalil P L Vq b ) 

+ B 3qaq b i% P L U ti ) («*. P^q b ) + B§ oft (u £j PrU^ ) (u qa P R Vg b ) 

+ ^ aqb (u ej P L u ei ){ug a P R v qb ) +Bl^ aqb (ue j P R ue i )(u qa P L v qb ) 

+ B iq a q b K ^^ P LU tl ) (U 9a a^P L V qb ) 
+ 0»vP R Ut % ) {Uq a a^P R V qb ) 

The very rich structure of the box-diagram amplitude is a consequence of the Fiertz 
transformation of the terms with product of lepton-quark and quark-lepton vector and 
axial-vector currents. All currents permitted by the Dirac algebra do appear. Each box- 
amplitude form factor has a chargino (C) and a neutralino (N) contribution. 

B L, R = B (N)L,R B (C)L,R , 2 2g v 

iq a qb iq a qb iq a qb ' v > 

= B^ ] n L ' R + B^ R . (2.29) 

Aqaqb Aqaqb oq a qb K ' 

Here and in the following equations q a qb assume the values appearing in the sum in Eq. 
(2.27). Neutralino contributions read 

B ^aqb = ^(^ A ^ fl ,<,<)iV2S*iVg^^ 

+ \d {M^,Ml v rnl x ^ (2.30) 

„{N)L _ !j (m 2 M 2 2 m 2 \ N R(e)* N R(e) N L(q> N L(q) 
B2q a q b ~ ~^\ M x\i M xV ex ' qY > iAX jBX bAY ^ aBY 

~ lMMl v M^,ml x ,ml)M^M^ B N^ (2.31) 

1 A7 -R(e)* ^(e) »rfl(?)* Ari(g) \ / 9 o 9 \ 

~~ 2 iAX ly jBX ly bAY ly aBYfi 

t5(N)L _ , , M 2 M 2 ^2 2 \f 1 AT R(e)* AT L(e) AT L(q)* ArR (q) 
B 3q a q b ~ a2(M-o A ,M xO B > m e x > m q Y ) \ ~ ^ N iAX N jBX N bBY aAY 

1 AT R(e)* AT L(e) AT L(q)* AT R ^\ (0 11\ 

~ 2 i AX Iy jBX ly bAY J W/' [4.66) 



AT R(e)* AT L(e) AT R(q)* A jL(q) 
ly iAX jBX bAY ly aBY f > 



} , (2-34) 



Cif = (< = l,-,4) (2.35) 



#3g a g b - B 3q a9b \L~R- (2.3b) 



(2.27) 
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The chargino contributions are 



„{C)L _ 1 , (M 2 M 2 2 ^2 s n R(e)* n R(e) n R{q)* n R(q) r 
°lqaq b ~ J 2 ^ Xl' Ya' i>x ' q' v > iAX L/ jBX°bBY a AY V 



4 v x A ' x B Ux i'y' 

+ ^(A^.A^^.m^^ V, (2-37) 

1 



2 v X A X b x 1y' *A x b 



= v 

- ^(M|_,M|_,<,^^ ^ (2 .3 8 ) 



1 „R.rf>.">* „r/ R 1 ^R/o 1 )* ^(nl ,. "I , n 



°3qaq b - yT ' y n ' x ' q'v \ O^iAX jBX bBY a AY °qd 



Xa 1 Xb 1 q'v' I 2 

R(e)* r ,L(e) n L(q)* n R(q) 



2°iAX ^jBX^bAY °aBy <V 



}, (2.40) 



r R{e)* r L(e) r R(q)* r L(q) <■ \ (r> a 1 \ 

~ °j.4X jBX bAY aBY °9"J ' l z - 4i J 

= (* = V", 4) (2.42) 

where for g = g' = d(n). The sum over paired indices is assumed. The Kronecker 

function 5 qu [5 qu ] denotes that q quark is one of up (u, c or t) quarks [one of the down 
quarks] . The loop functions do and d>2 are evaluated by neglecting the momenta of incoming 
and outgoing particles. They are listed in Appendix D. 



3. Amplitudes and branching ratios 
3.1 Hadronization of currents 

The effective Lagrangians (the matrix elements) for photon and Z-boson part of the am- 
plitudes for SL LFV lepton decays comprise vector and axial-vector currents, but box 
amplitude contains all possible quark currents permitted by Dirac algebra, that is scalar-, 
pseudoscalar-, vector-, axial-vector- and tensor-quark currents. To perform calculation of 
charged-lepton SL LFV decays rates these currents have be converted into meson currents 
comprising mesons that appear in the possible final products of the charged-lepton SL 
LFV decays we study here. The hadronization procedure we use here is not exact in the 
sense that we do not include the sea-quark and gluon content of the meson fields, but 
it is precise enough to give much better than order of magnitude decay rates of the pro- 
cesses considered here. The quark content of the meson states is given in Appendix E. The 
hadronization of axial- vector current is achieved through PCAC (see e.g. ([j34[, |35| , [3t|); 
for normalization of pseudoscalar coupling constants used here and for further details see 
@ ) . The hadronization of the vector current is achieved using vector- meson-dominance as- 
sumption (see ( p5| , |36|); for normalization of the vector meson-decay constant and details 
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see |37| , 0, || ) . Hadronization of scalar currents is achieved by comparing the quark sector 
of the SM Lagrangian and the corresponding effective meson Lagrangian |38|] (for applica- 
tions in content of LFV and details see ||]). Hadronization of the pseudoscalar current is 
obtained by the same procedure as for the scalar current. The results obtained by using 
this procedure is equal to the result obtained by using equation of motion for current-quark 
masses (see e.g. |39|) and results for hadronization of axial- vector current up to the differ- 
ence of up and down quark masses or up to the difference of pseudoscalar decay constants. 
The hadronization of the tensor-quark currents is obtained by comparing the derivative 
of tensor-quark current with vector-quark current and using the equations of motion for 
current quark masses. The difference between terms, one containing the derivative of the 
incoming quark field and the other containing derivative of the outgoing quark field, have 
been neglected. The error expected from this approximation is proportional to the amount 
of breaking of 5 , C/(3)fl avour symmetry. The tensor currents are proportional to the cur- 
rent quark masses, and therefore give smaller contribution than the other quark currents. 
Therefore, the error introduced by this approximation in the total SL LFV amplitude is 
small. 

Here we summarize the basic quantities needed to describe the hadronization of quark 
currents. 

1. The pseudoscalar meson decay constants [^3| (fp, P = 7T°, 77, rf, K°,K°); 

2. The constants 7^ |57) (V = p°, <ft, uj, K*°,K*°) defining the vector meson decay 
constants (/y ~ rriy/jv)] 

3. The mixing angles Op and Qy E3| defining the physical meson- nonet states in terms of 
the unphysical singlet and octet meson states; 

4. The parameter r |38j (m u , m,d and m s are current quark masses), 

2ml+ _ 2m 2 K0 _ 2m 2 K+ 
m u + m d m d + m s m u + m s ' 

that appears in the hadronization procedure for scalar and pseudoscalar currents. 

Having the identification of the quark currents with the corresponding meson cur- 
rents, achieved by the above hadronization procedure, one can write down the effective 
Lagrangian as a sum of terms with an incoming lepton field li and outgoing lepton field £j 
and pseudoscalar meson (P) or vector meson (V) field(s). This Lagrangian directly gives 
the amplitudes for the li — > ljP(V) processes. Amplitudes with a pseudoscalar meson in 
the final state contributions come from the pseudoscalar and axial-vector coupling part of 
the effective Lagrangian, while the amplitudes with vector mesons have vector and tensor 
coupling contributions. Only the scalar coupling gives no contribution to the one- meson 
processes in the final state, li — » £jP(V). They contribute only to the processes with two 
pseudoscalar mesons in the final state, ii — > ijP\P2- 

3.2 Vector-meson pseudoscalar-meson interactions 

The processes with two pseudoscalar mesons in the final state are generated by the scalar- 
quark-current part of the effective Lagrangian, and vector-quark- and tensor-quark-current 
of the effective Lagrangian. The scalar-quark-current part of the effective Lagrangian 



produces two pseudoscalar fields directly. The vector-quark- and tensor-quark-current parts 
produce a resonant vector meson state (V), which decays into two pseudoscalar mesons 
(P). The VPP interactions necessary for description of the VPP interactions appearing 
in the charged-lepton SL LFV decays are described by the part of the meson Lagrangian 
containing these VPP vertices || 

C VPP = -^{^(27^ vr- + K+d»K- - K°d,K°) 

+ K°*^(V2TT-d^K + - ir°d p K° - V3c P K°d^V ~ V3s P K°d^')} + ■ ■ • -(3.2) 

This Lagrangian is a part of the nonlinear (U(3)l x U(3)p)/U(3)y symmetric sigma-model 
Lagrangian. U(3)y symmetry corresponds to the vector mesons in the linear realization of 



the gauge equivalent (U(3)l x U(3)p) g \ ^ a _i x U(3)y linear sigma-model f4Q, 41]. One can 



include the (U(3)l x U(3)p)/U(3)y breaking terms, too [41]. That was applied to SL LFV 
tau-lepton decays in Ref. ||, but for the estimates of the charged-lepton SL LFV decays 
it is unnecessary complication, and we will not consider it here. 



From Eq. (3.2) one can read of the cyp 1 p 2 couplings in terms of g pn7T coupling. For 
instance, c p o K+K - = ^g p7nr . 

When the amplitudes with vector- meson resonance (s) are formed, the square of the 
vector meson mass in the rriy/'yy, appearing in every vector meson decay constant, has to 
be replaced with (m v — imyTy) /^/y , where Yy is the decay width of the vector meson (8|. 

3.3 Charged-lepton SL LFV with one meson in the final state 

Now we can write down all amplitudes we are interested in. (Notice that the axial-vector 
mesons (^4) are not included. They decay into three pseudoscalar mesons, and therefore the 
amplitudes are much more complicated. New vertices with A — V — P couplings should be 
included, kinematics is much more involved etc. The scalar mesons are also not included. 
They also lead to complications.). They are given as follows. 



+ q 2 1 2 7 L + g 2 1 2 7 R 



ijV 
IR 



+ ia^Pm^il + w^Pm^l 



utietf. (3.3) 



^ + B^)+j,Pr(P^+B 



+ 



P L B^ + P R B^ 



(3.4) 
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The form factors in Eqs. (|3.3j ) and (^J) are defined as follows: 



P. 



ijV 
07 L,R 



(« = 1,2), (V = p°,4>,u) 



V 



ijV 
ZL,R 



9 2 Z V L,R m y k V 



,2 J2 ' Z 



V2~fv 



Z ! 



ijV _ m v 1 



(V = p°,<t>,u,K*°,K^), 



a 2L,R 



-2V2 

7y 



P 



ijP 



+ k gd (m s - m d )B^: d 



T L z > R {-V2f P )k p z , 



m%c 2 



(V = K*°,K*°), 

(P = 7T ,7/,»/) 



&ilr = 8 L ,nk/2fp) 



~B\ku + B 2Uu) + ^dd 



Ids *' 2ds 



Bf B = ^(~)(v^/r^ 



ldd + 2dd 



■2L.H " S £,M- 2")(V2/ P J 2 



(P = Tr ,7], 7]', K",K°) 
P 



„L,R J5L,R\ , 



°3dd + °3dd 



(P = ^°,r ? ,r ? ',K ,^). 



(3.5) 
(3.6) 



(3.7) 

(3.8) 
(3.9) 



(3.10) 



(3.11) 



In Eqs. ( 3.1 0| ) and (3.11) sl = 1 and sp = —1. The constants , kg, k%, and fe^^ are 
defined in Appendix E. 

A branching ratio for the processes — > £j-P with unpolarized initial and final particles 
reads 



B(4 -> ^P) 



111 A2 (m? , mj , rrip) 
8tt mf IY 2nii 



I 7 ZL + °1L I + I ' ZR + °LR I 



+ [\*£\ 2 + \*#£\ 2 )ip2 



(TTl + Bi j L P )(V^ + Bgf )• + c.c.) m jmi i P3 
+ {(P$Z + BgW)* + (7^ + B^)(B^T + c.c.) myi« 
+ ((PSl + BgW)* + {vTr + B^){B^T + c.c) m iiP5 



- 13 - 



+ 



((BV L P )(B^y + c.c.)m imj 



where I\ is a total decay rate of the lepton k and 

1 



i pi 



IP2 



((m 2 — rn 2 ) 2 — (mf + rnj)mp) 



(mf + m 2 — mp), 



i P3 = m p 
ip4 



/ 2 i 2 

(m i + m P — m 



2 V "' 



/ 2 2 2\ 

(mj — m P — rrij), 



and 



X(x,y,z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 



(3.12) 



(3.13) 



(3.14) 



A branching ratio for the processes — > ^-V with unpolarized initial and final particles 
reads. 



111 A2 (mf,mj,m^) 
87r mf 2mj 

(I^l + ^ + b il\ 2 + \Kr + nl + b ir\ 2 ) ivi 



+ 



V, 



ijV 

27L 



rn 



2 + °2L 

y 



+ 



V. 



ijV 
2-yR 



rn 



2 ^ °2_R 

y 



«y2 



+ 



+ 



p. 



ij-y \ 

2 7 L + ^ 



V 



m 



2 ^ °2_R 



+ c.c. (—myrriirrij) 



v 



ijVy. 



v. 



ijV 

27L 



m 



y y 



p. 



ij-y 

2 7 i? 



+ #2iT ) + cc. I {mjivs) 



m 



V 



V. 



ijV 
2fR 



1L 



rn 



2 + B 2R 
V , 



+ I' 1-yi? ^ ' ZR + °1 



P. 



y-y 



m 



+ c.c. (miZy 4 ) 



V 



(3.15) 



where 



iyi = [ m y( m ? + m i) + ( m i ~ m j) 2 ~ 2m y] > 



- 14 - 



W2 



, o 2\2 ^~ 2 / 2 2\ ^ 4 

(m< -m^) - -m y (m i -hm,-) - -my, 



1 



«V3 = 2 l m i 

1/ 2 
«V4 = ^ ( m i 



m 2 + my) 



2 2 \ 
rrij — m v ) 



(3.16) 



3.4 SL LFV decays of a lepton with two pseudoscalar mesons in the final state 

Amplitude for a general £j — ► ljP\Pi decay rate is a sum of scalar-coupling contribution 
and resonance contributions (coming from vector- and tensor-coupling contributions) 

,PiP 2 .• A Ji-*£jPlPi , ■ A A ti-*liP%P% 



where 



M, 



iM l 



mi. 



iM%7 £jPlP2 + iM l { 



(3.17) 



711 o 



H)Pl 



9 

m| 



H)Pr 



+ E^Ha^PLipt -p 2 Yq v + E p R p Ha ilu P R ( Vl 
= iu tj [P L Af PiP2 + P R A? PiP2 ] u k . 



(3.18) 
(3.19) 



In Eq. ( |3.18 ) D P ^ P R and E P ^ P p are form factors built from the trilinear cvp 1 p 2 couplings 
(defined by the Lagrangian (|3.2| )), normalized vector-meson propagators, 



my - imyTy 
7 2 — my + imyTy 



(3.20) 



and form factors for 
meson, e.g. 



jV processes divided by the the mass of the resonant vector 



V 



ijV 



V 



ijV 



m 



V 



The expressions for the D p ^ p p and E^j^ p form factors are 



P1P2 



D PiP 2 



E PiP 2 
^IL,R 



l7i,R ' ZL,R °1L,R 



111 



v 



imyTy 



E 

V 



n 



ijV \ 
2 7 L,R 

2 °2L.ii 



mf 



q 2 — m y + imyTy 
- imyTy 



CVP1P2, 



f — m v + imyTy 



Cy Pl p 2 . 



(3.21) 

(3.22) 
(3.23) 



The sum goes over neutral vector mesons only (V = p°, (j), u, K*°, K*°). The coefficients of 
the non-resonant part of the amplitude, the constants A\ Pl p 2 , are defined as the coefficients 
of the P1P2 product of fields contained in the matrix- valued operator 



1 E 



(IT 



kbqa(^3q aqb + ^3q a q b ) 



(3.24) 
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where II is a matrix of pseudoscalar fields 



n 



\/2vr- 



V / 2vr+ 
^ + 7=3™ + ^ 



^0 _^_^V2, 



(3.25) 



For example, 



I f B L,R , T5L,R L)R -pL,R 



(3.26) 



Having the amplitudes one can easily evaluate the branching fractions. We assume that 
incoming and outgoing particles are not polarised. 



B(£i 



1 1 



I 



(rrii-mj) 2 



ds 



(2vr) 2 Y U 32m? „ 

'(mi+m 2 ) 2 
I I Ami-mj) 2 

dsn 



12 



rn a :c 
\#2 



mi n 
J2 



ds j2 \M 



e i ^e J p 1 p 2 \2 



(2vr) 2 32mf 



i J(mi+m 2 ) 2 



<i P2 i 2 + i^ P2 i 2 Ui 



+ ii^i 2 +i^n/2 



'IP 



pPlP 2 |2 



+ (\A 



IPiP 2 \ 



R 



V* 



+ 
+ 



+ {(A L lPlP2 ){A? PiP2 y+c.c)h 



+ 
+ 
+ 
+ 



P\P2\* 



+ C.C.J I 7 

+ cc. ) 7 8 



pfi P2 )(^ P2 )* + p p A P2 )(^r 2 : 

P P i P2 )(^ P A P2 )* + P P A P2 )(^ P2 ^ 

«2 P2 )(Af PlP2 )* + (D^)(A? PlP2 y + cc.) i 9 
(D? L P2 )(A? PlP2 Y + «/ 2 )(Af PiP2 )* + cc.) I 10 



+ ((^ P2 )(^p iP2 )* + (^T 2 )(Af PlP2 )* + c.c.)h 



PIPn 



IB 



The Mandelstam variables are defined as s a b = (p a — Pb) 2 , e.g. S12 = (pi — P2 
kinematical bou 
I integrals read, 



kinematical bounds on the Mandelstam variables s 1 J^ n and s^ ax are well known 



(3.27) 

. The 
The 



h = 2s 2 2 + 2s] 2 ~(e 1 + e 2 ) + (2eie 2 - e 3 e 4 )l, 

J 2 = _2s 2 2 (ei ) + 2sj2(e 5 ei + egeio - e 6 e 7 - e 8 e 7 ) 

+ 2((e 5 e 6 e 7 + e 8 e 9 e 7 - e 5 e 9 ei - e 8 e 6 en) + e 3 (eneio - e 7 ))T, 



- 16 - 



^3 = e 3 l, 

I 5 = rmmj (ei eu - e 2 )T, 



-^6 = mim jl 



I 7 = m,-(-ei 2 e 6 )l, 
J 8 = mi(ei 2 e 8 )T, 



(sj2 + (e 2 )l) , 



-Z10 = «^ (sj 2 + (ei)l) , 
111 = -Sj2(e8 - e 6 ) + (e 8 e 9 



where 



„n 

s i2 



. ma x 
V2 



J2 



The quantities ej read 



ei 


= -es - 


Ti2 — mi 

e 8 

S12 


e 2 


= -eg - 


m 2 — mi 

e 6 

S12 


e3 


1/ 2 
2 1 


+ rrij- s 12 ) 


e 4 


= en - 


(ml — m 2 ) 2 


S12 


es 


= m 2 + 


f 1 2 

-(m^ + ro j - 


ee 


= 2 K 


- m 2 + s 12 ) 


e7 


= m 2 — 




es 


1/ 2 
= 2 K 


- m 2 - S12) 


eg 


= mf + 


I/O 

-(mi +m^ - 


eio 


= S12, 




en 


9 9 

= 2m x + 2m 2 — si 2 , 


ei2 


= -e 4 . 





4. Minimal SO (10) model and its predictions 



(3.28) 
(3.29) 



(3.30) 



Now we are ready to estimate the SL LFV processes. In order to perform a concrete 
evaluation for the SL LFV processes, we need an information on the Yukawa couplings. 
In this paper, we make use of the minimal 50(10) model, as an example which gives a 
precise information for the neutrino-Dirac- Yukawa couplings. We begin with an overreview 



of the minimal SUSY S'O(IO) model proposed in [22] and recently analysed in detail in Ref 



1 23, 24, 25, 26, 27 1. Even when we concentrate our discussion on the issue how to reproduce 
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the realistic fermion mass matrices in the SO (10) model, there are lots of possibilities for 
introduction of Higgs multiplets. The minimal supersymmetric SO(W) model is the one 
where only one 10 and one 126 Higgs multiplet have Yukawa couplings (superpotential) 
with 16 matter multiplets. Therefore, the quark and lepton mass matrices can be described 
as 

M u = C W M W + Ci2 6 Mi2 6 , 

M d = M 10 + M 126 , 
M D = c 10 M 10 - 3c 12 6Afi2 6 , 
M e = Mio -3Mi2 6 , 

M R = c R M 126 , (4.1) 

where M u , M d , Md, M e and Mr denote up-type quark, down-type quark, neutrino Dirac, 
charged-lepton and right-handed neutrino Majorana mass matrices, respectively. Note that 
all the quark and lepton mass matrices are characterized by only two basic mass matrices, 
Mio and Mi 2 6> and three complex coefficients cio, C126 and cr. 

The mass matrix formulas in Eq. ( |4.1| ) lead to the GUT relation among the quark and 
lepton mass matrices, 

M e = c d (M d + kM u ) , (4.2) 

where 

3cio + C126 / . „n 

Cd = , (4.3) 

ClO — Cl26 

k = . (4.4) 

3ci + C126 

Without loss of generality, we can start with the basis where M u is real and diagonal, 
M u = D u . Since Md is the symmetric matrix, it can be described as Md = V£ KM Dd Vq KM 
by using the CKM matrix Vckm and the real diagonal mass matrix Dd- 1 Considering 
the basis-independent quantities, tr[MjM e ], tr[(MjM e ) 2 ] and det[M=M e ], and eliminating 
\cd\, we obtain two independent equations, 

f tr[M e f M e ] ^ 2 ti[(M e ] M e f] 

(4.5) 



nij + m 2 + rn^. J m\ + + ' 

(4.6) 



/ tr[M e f M e ] \ = det [M e f M e ] 
I ml + m 2 + m 2 J m\ m 2 m 2 ' 

where M e = V^ KM V^ KM + kD u . With input data of six quark masses, three angles and 
one CP-phase in the CKM matrix and three charged-lepton masses, we can solve the above 



general, M d = U* D d [/ f by using a general unitary matrix U = e ia e ipTa e ilTs VcKMe i ' 3 ' T3 e i ' y ' Ts . We 
omit the diagonal phases to keep the number of the free parameters in the model as small as possible. 
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equations and determine k and \cd\, but one parameter, the phase of c^, is left undetermined 



1 23]. The original basic mass matrices, M\q and M\26, are described by 



^ Z + \^L v$KuDd yi Ku + \«\f± Dnt ( 4. 7) 



M 126 = 1 l l dlela V6KM Dd 4m " l -^^D u , (4.8) 

as the functions of a, the phase of cy, with the solutions, |c^| and k, determined by the 
GUT relation. 

Now let us solve the GUT relation and determine |c^| and k. Since the GUT relation of 
Eq. ( |4.2[ ) is valid only at the GUT scale, we first evolve by renormalization group equations 
(RGE's) the data from the weak scale to the corresponding quantities at the GUT scale 
for a given tan (3 Then we and use them as input data at the GUT scale. Note that it 
is non-trivial to find a solution of the GUT relation, since the number of free parameters 
(fourteen) is almost the same as the number of inputs (thirteen). The solution of the GUT 
relation exists only if we take appropriate input parameters. By using the experimental 



data at the Mz scale 42], we get the following absolute values for charged fermion masses 
(in units of GeV) and the CKM matrix at the GUT scale, Mq, with tan (3 = 45, m s = 0.072 
and 5 = 1.518: 

m u = 0.00103 , m c = 0.299 , m t = 133 
m d = 0.00170 , m s = 0.0263 , m b = 1.55 



0.000411 , rrif, = 0.0868 , m T = 1.69 



and 



Uckm 



/ 0.975 0.222 0.000146 - 0.00279z 

-0.222 - 0.00012H 0.974 + 0.000129i 0.0320 
y 0.00697 - 0.00272i -0.0312 - 0.000626i 0.999 



in the standard parameterization. The phases of the fermion masses are not determined 
by the diadonalization procedure. Here the masses are chosen be real. The signs of the 
fermion masses have been chosen to be equal, — for m u , m c , rrid and m s , and + for mt 
and rrif,. So determined masses at the GUT scale are used as input parameters in order to 
solve Eqs. (|4.5| ) and (^^). As an example, here we show one of two solutions, 

k = 0.0134 - 0.000791 i , 
\c d \ = 6.39. (4.9) 

Once the parameters, \cd\ and k, are determined, we can describe all the fermion mass 
matrices as a functions of a from the mass matrix formulas of Eqs. ( |4.l[) , ( [4.7|) and (Un- 
interestingly, in the minimal S'O(IO) model even light Majorana neutrino mass matrix, 
Mjy, can be determined as a function of the phase a and cr through the seesaw mechanism 
M v = -M^M^Md. 
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Now we give an example of the neutrino-Dirac- Yukawa coupling matrix which can fit 
the GUT relation. In the basis where both the charged-lepton and right-handed Majorana 
neutrino mass matrix are diagonal with real and positive eigenvalues, the neutrino Dirac 
Yukawa coupling matrix at the GUT scale for fixed a = 3.223 [rad] is found to be 



Y v 



( 0.000310 + 0.00348i -0.000894 - 0.000249i 0.0447 + 0.053H ' 

0.00590 - 0.0103i -0.0164 - 0.0427i 0.308 + 0.116i | . ( 1.10) 
^ 0.00215 + 0.00126i 0.0558 - 0.0559i -0.381 + 0.604i 



Using these Yukawa coupling constant matrices, we proceed with numerical calculations. In 
evaluating the LFV branching ratios, we first solve the RGE's for the soft SUSY breaking 
parameters and the Yukawa couplings in the MSSM to determine the masses and mixings 
for the SUSY particles. Then we input these data into the formulae presented in the 
previous sections. 

In the following, we list the input parameters what we used for the hadronization 
processes. For the pseudoscalar meson decay constants, we take as an input the following 
values [43 1, 



Uo = 0.119 [GeV], /„ = 0.131 [GeV], = 0.118 [GeV], (4.11) 

and for the vector-meson decay constants we use the values, extracted from the vector 
meson decay, V — ► e + e~, 

7p o = 2.518, 70 = 2.993, 7w = 3.116. (4.12) 

The mixing angles between singlet and octet states for the vector mesons and for the 
pseudoscalar mesons that we use are [fi3|l , 

V = 35°, 6 P = -17.3°. (4.13) 

In order to investigate the dependence of SL LFV branching ratios on model parameters 
we plot their dependence on m^: r — > en , r — ► erj, r — > erf , r — > /U7r°, t — ► fir], and 
r — ► \irf graphs in Fig. 1, r — ► ep°, r — > fj,p r — > e<f>, r —>■ \i<$> r — > eu>, and r — > iiu graphs 
in Fig 2, r — » eir + ir~ , r — > fj,ir~ir + , r — > eK°K , and r -> iiK°lC , r -» eK + K~ , and 
r — > iiK + K~ , graphs in Fig 3, and r — > e 7 and r — > //7 graphs in Fig 4. In Fig 5. we plot 
the TOf fi dependence of the quantities used to determine the alowed region of mf R mass, 
the MSSM contribution to the anomalous magnetic moment Aa^ and \x — > e 7 branching 
ratio, The parameters are chosen so as to satisfy the WMAP constraint on the cold dark 



matter (CDM) relic density [45|, 



n CD M^ = 0.1126. (4.14) 
We can transmute this value into the approximate relation between mo and M x / 2 , such as 



m [GeV] = ^M 1/2 [GeV] + 150 [GeV], (4.15) 



- 20 - 



BR(x 




-14 

1.8x10 












-14 

1.6x10 


















1.4x10 14 


















-14 

1.2x10 


















1.0x10 14 


















0.8x10 15 



















200 250 300 350 400 
m- [GeV] 



BR(x [L%°) 

2.5xl0" 12 



2.0x10 



1.5x10 



1.0x10 




200 250 300 350 400 
m- [GeV] 



BR(x^ eq) 

3.0xl0" 15 



2.5x10" 



2.0x10 



1.5x10 



15 



15 



200 250 300 350 400 
m- [GeV] 



BR(x -> |xti) 

13 



4.5x10 



3.5x10" 



2.5x10" 



1.5x10" 




200 250 300 350 400 
m- [GeV] 



BR(x -> eq') 

3.0xl0" 15 



2.5x10 



2.0x10 



1.5x10" 




200 250 300 350 400 
m- [GeV] 



BR(x -> [IT]') 

4.0x10"' ' 



3.5x10" 



3.0x10 



2.5x10 



13 



13 



2.0x10 



1.5x10 



200 250 300 350 400 
m- [GeV] 



Figure 1: The branching ratios for SL LFV decays r — ► en / erj / er}' / /j,tt° / /ir] / \ir( as a function of 
mass of lightest charged sfermion mf R . As an input we have taken tan (i = 45, n < 0. and A = 0. 



for tan/3 = 45, and Aq = 0. The maximal values for the branching ratios of the r — > e//x7r' 
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Figure 2: The branching ratios for SL LFV decays t — > ep° /e<p/eoj/ up / fi(f>/ fiuj as a function of 
mass of lightest charged sfermion mf a . The input parameters are as in Fig. 1. 



processes are found to be 

BR(r -► evr ) ~ 1.7 x 10~ 14 , (4.16) 
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Figure 3: The branching ratios for SL LFV decays r — ► 

en + n~ /eK°K /eK + K~ /^n + -K^ / [iK + K~ /^iK°K as a function of mass of lightest charged 
sfermion m? B . The input parameters are as in Fig. 1. 



BR(t ^vr ) ~ 2.4 x 10 



-12 



(4.17) 
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Figure 4: The branching ratios for decays r — > ej/fij as a function of mass of lightest charged 
sfermion rrif R . The input parameters are as in Fig. 1. 




Figure 5: The anomalous magnetic moment and branching ratios for /i — > e7 as a function of mass 
of lightest charged sfermion mf a . These quantities restrict the region of rrif R values. The present 
experimental upper limits are represented by horisontal red lines. The input parameters are as in 
Fig. 1. 



and for the r — > e//j,rj processes BR(r —>■ e/iiif) ~ 0.15 x BR(r — ► e/[nr°), with suitably 
chosen CMSSM parameters which can realize the neutralino dark matter scenario by the 
WMAP data. It can be realized with the following set of parameters: tan (3 = 45, \x > 
and Aq = 0, with Mi/ 2 = 600 [GeV], tuq = 343 [GeV]. This parameter set can also predict 
the muon g — 2 within the range of the recent result of Brookhaven E821 experiment and 
also provides the r — * ^7 and — ► erf branching ratios close to the current experimental 
bound H . The ratio between two processes r — * e/fiir and r— > e/ fir] is & result of the 
dominance of the Z-boson-penguin amplitude in these processes, and reflects the difference 
in the form factors and the mixings between the singlet state and the octet state of the r\ 
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x ^ + ^ ~ 0.15. (4.18) 



mesons, 

BR(r^e//xvr) V/^ " VV3 ' \/6, 

We can also see the correlation between the branching ratios for the processes r — > e/p p° 
and r — > e//x7 as BR(r — >• e/pp°) ~ 5.5 x 10~ 3 x BR(r — > e/fij). It can be estimated on 
the basis of the photon-penguin-amplitude dominance in these amplitudes, giving 



BR(T^e/,. 7 ) 2 Vv 

When we impose the constraints from discrepancy of r and e + e~ data in the muon g — 2 
measurements [46] and from upper limits on the p — > e^f branching ratio, we obtain that 
the model permits the of rrif R values satisfying, 

m fR > 204 [GeV] . (4.20) 

The lower bound comes from the muon g — 2 constraint. The p — > cy gives also the lower 
bound but it is below the lower bound from the muon g — 2 constraint. The g — 2 curve 
has uprising behaviour above m? R = 340 [GeV], but at mf R = 1000 [GeV] it is almost 
independent on m? R and has a value 5.4 x 10~ 10 , slightly below the present experimental 
g— 2 uncertainty. Therefore, one can expect that the improvement of the g— 2 measurements 
will give the upper limit on nif R , too. Using the lower bound on rrif R values one can find 
the theoretical upper bounds for all leptonic and SL LFV branching ratios. Leptonic and 
dominant SL LFV deacays are given in Table I. 

5. Summary 

The evidence for the neutrino masses and flavour mixings implies the non-conservation of 
the lepton- flavour symmetry. Thus, the LFV processes in the charged- lepton sector are 
expected. In supersymmetric model based on the minimal S'O(IO) model, the values for 
the rates of the LFV processes are generally still several orders of magnitudes below the 
accessible current experimental bounds. In this paper, we have presented the detailed 
theoretical description for the SL LFV decays of the charged leptons with one or two 
pseudoscalar mesons or one vector meson in the final state. Also, some previous formulae 
have been corrected. The 7-penguin amplitude is corrected to assure the gauge invariance, 
the Z-penguin amplitude is corrected, new box contributions to the box amplitude has 
been found and previously neglected terms are given. 

To evaluate the decay rates of the LFV processes within the MSSM, the parameters 



and the LFV interactions of the MSSM have to be specified. It has been shown [23] that 
the minimal SUSY 50(10) model can simultaneously accommodate all the observed quark- 
lepton mass matrix data involving the neutrino oscillation data with appropriately fixed 
free parameters. In this model, the neutrino-Dirac- Yukawa coupling matrix are completely 
determined, and its off-diagonal components are the primary source of the lepton-flavour 
violation in the basis where the charged-lepton and the right-handed neutrino mass matrices 
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Theor. upper bound 


Exp. upper bound 
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— ► e7 
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1(T 13 


1.2 x IO -11 
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— ► e7 
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1(T 12 


1.0 x IO -9 
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-» ^7 


4.1 x 


10 -io 


4.5 x IO -9 




T 




1.7 x 


10 -14 


1.9 x IO -7 




T 




2.4 x 


io~ 12 


4.3 x 10~ 7 




T 


— > ery 


2.9 x 


io- 15 


2.3 x 10~ 7 




T 


-* W 


4.2 x 


io- 13 


2.3 x IO -7 




T 


— ► er/ 


2.8 x 


io- 15 


10 x 10~ 7 




T 


-> /if/ 


3.7 x 


io~ 13 


4.1 x IO -7 




T 




2.3 x 


IO" 14 


2.0 x 10~ 6 




T 




2.8 x 


1(T 12 


6.3 x 10~ 6 




T 


— > e0 


1.4 x 


1Q -14 


6.9 x IO -6 




T 




1.7 x 


10 x/ 


7.0 x 10 




T 


— ► ew 


1.6 x 


10 _ 










J. . O A 






T 






4.2 x 


1Q -14 


8.7 x 10~ 7 


T 






5.8 x 


io- 12 


2.8 x 10~ 7 


T 




eK°K° 


4.3 x 


io- 15 


2.2 x 10~ 6 


T 




HK°K° 


5.3 x 


io- 13 


3.4 x 10~ 6 


T 




eK+K- 


6.7 x 


io- 15 


3.0 x 10~ 7 


T 






5.1 x 


IO" 13 


11.7 x IO" 7 



Table 1: Theoretical upper bounds I — > £'j processes and dominant SL LFV processes. The upper 
bound is obtained from the muon g — 2 constraint. We referred the experimental data mainly from 
@ and partly from |43| . 

are real and diagonal. Using this Yukawa coupling matrix, we have calculated the rate of 
the LFV processes assuming the mSUGRA scenario. The analytical formulae of various 
SL LFV processes, are given. Using these formulae, 

we have numerically evaluated £j — ► £jP, ii — ► £jV and £j — ► ijP\P2 branching ratios. 
Among these, the branching ratios of r — > n r ) and /z — ► ej may be interesting for the near 
future experiments. The typical CMSSM parameters used in calculations can realize the 
neutralino dark matter scenario by the WMAP data. 
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A. Notation for the MSSM Lagrangian 

Here we summarize our notation necessary for defining the masses of the sparticles in the 
MSSM Lagrangian. 

V EE ^(^0)2 + {H 0f = lUA [GeV]; (A-1) 

and 

Then the charged fermion mass matrices are given by 

M$ = -l>sin/3, (A.3) 
M$ = Yi j vcoap, (A.4) 
Ml j = Y* j vcosp. (A.5) 

The mass matrix of the charginos is written as: 



£ = -(w r ,h u 



R 



H 



L +h.c, 



dL 



M 2 V2M W cos (3 
V2M W sin [3 \i 



(A.6) 



The mass matrix of the neutralinos is written as: 



C = -\[B U Wl,H° L ,H° L ) M^o 
( 



wl 



+ h.c, 



Mi -M z sm 9 w cos/3 M z sin 6w sin (3 \ 

M2 Mz cos 0w cos (3 —Mz cos 9\y sin (3 

Mz sin 9yy cos (3 M z cos 9w cos (3 —fi 

\ Mz sin 9w sin (3 — M z cos 9w svo. f3 — ^ o J 

(A.7) 



The mass matrices of the squarks are written as follows: 
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m 7l 



2 _ " l q 



m 2 - + MlM u -Alv sin - (i* cot /3 



" " » -A u v sin /3 - M u fjL cot /3 m| + M u i\4 - Mf cos 2/3 sin 2 W 

M 2 cos2/3(±- |sin 2 ^)l3x 3 \ 

§M§ cos 2/3 sin 2 6 w l 3x3 J ' 



+ 



2 _ / m| + M]M d cos /3 - M] /x* tan /3 \ 

mj ~~ y cos - M d /j, tan (3 m 2 - + M d Mj - Mf cos 2/3 sin 2 W y 

/M 2 cos2/3(-i + isin 2 ^)l3x 3 \ 

V -±M 2 cos 2/? sin 2 9 w l 3x3 J ' 1 

The mass matrices of the sleptons are written as follows: 

£ = -(mlr fyi - (mfr efe, 
m? = m| + -Mf cos 2/31 3X 3, 



mf 



' m| + M e f M e cos /3 - m] fi* tan /3 \ 

^ A e v cos /3 - M e ix tan /3 m| + M e i\l\ - Mf cos 2/3 sin 2 6 W J 

/M 2 cos2/3(-i + sin 2 ^)l 3 x3 \ 

y —Mf cos 2/3 sin 2 6*^/13x3 y ' 1 

They are diagonalized with unitary matrices as follows: 

OrM^ 0\ = diag(M~-,M~-), 

0* N M^o O j N = diag(M~o , M~o , M~o , M~o ) , 

U f m j U f = diag(m m~), (/ = u,d,e), 

^ ™~ Ul = diag(m| 1 , m\ , m\ ) . (A. 10) 



B. Lagrangian for fermion-sfermion-gaugino/Higgsino interactions in MSSM 



The LFV interactions in the MSSM include the fermion-sfermion-gaugino/Higgsino ver- 
tices. These vertices, and corresponding coupling constants (Cf^x^ , N^J^^ , / = v, e, u, d) 
are defined by the following Lagrangian 



C = Ui 



R{u) 



iAX r R 



X + A dx + d t C^Pl + CV^Pr 



R(d) 



X A U X 



+ u i C iAX )p RXA e X + di 



R(e) 



iAX 



AX J 



+ Ui 



N-a^Pl + NITx'Pr flfix + di K^Pl + N^x'Pr X\dx 



T R(u) 



L(d) 



T R(d) . 



+ ws$p«&vx + * [n^pl + ^s^; 

+ h.c. 



Xa^x 



= UiP L x\dx 



m. 



v^Mwsin^ 



(0 R ) A 2m)xi 
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rn 



\/2¥ w cos (3 
g{(0 R ) A1 (U*,) Xi } 
9 f m u 



\/2M w cos (3 



+ u l P R x\d x g\(0 L )Ai(U})xi + ™ dk S V CKM ) kj m KM ) ki {0 L ) A2 {Ul) X j +z 
L I V2M w cos/3 a 

+ 3ii^xA^ \a{ ^™ dl R {0l) A2 {ui) x \ (VcKuhj 

I l v 2Mw cos p J 
+ ^P«xI% 9\{0* R )AiW)xj ~ " J . - (0* R ) A2 (Ut) x , j+ z (Vg KM ) 

+ ^xlex[ 5 ((OLU(t/|)x j + q (Ql)^(^|)^-+3 j (E^Ns)ij 

+ eiP L x A 9 x 

+ eiP R x A i>x 

+ UiP L x A ux 
+ UiP R x° A u X 
+ diP L x\d x 
+ diP R x\d x 
+ 



7f{i^i^ (0jv)A4( ^ k ' i+3 + 

5 f "id 



(CV)A2 + - tan ^(Oat)ai 



g / 



v/2\ M w cos/3 



(OAr)A3(^|)x,j+3 



(CVU2 + -tan6»vi/(OAr)Ai 



(Vckm) 



+ ^-PrxVx ^= [(Oat)a2 - tw-0 w {O N ) A1 ] (U~) XJ (U M Ns)ij 



+ eiP L x A e X ^{- M ^ e c l s(3 (°*N)A3(U!)xi + 2t aa w (O* N ) A1 (Ul)x, i+ 3j 



+ e t P R x° A ex -^{- A/r mei R (0 N ) A3 (U!) x , i+3 + [-{0 N ) A2 -t W 8 w {0 N ) Al ](Ul) Xl 
V 2 ^ Mw cos p 

+ /i.e. (B.l 

C. Trilinear interactions of fermions or bosons with Z°-boson or photon 

The interactions of Z°-boson or photon with any fermion or any boson follow from SU(2) l X 
U(1)y gauge symmetry. 



(C.l) 



(C.2) 



The interaction Lagrangians are flavour-diagonal in the weak basis. In the mass basis the 
interactions with photon remain diagonal, because all mixed states A must have the same 
charge, Q A . On the other side the interactions with Z-boson, which depend on charge 



* = -»£ (»* & a) \z> (/f-L - Q f>-k) + (w) 
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(Qa) an< i third component of the weak isospin (/3A) are not in general flavour-diagonal in 
the mass basis, because the mixed states may have different I^a values. 

For LFV processes the interaction of photon and Z-boson with charginos, neutralinos 
and sfermion fields is needed. The Lagrangians for the corresponding weak-basis fields is 
easily written knowing the charges and J3-S of these fields. After transformation from the 
weak-basis to the mass-basis the following interaction Lagrangians are obtained 



C 



(E L A { pP L + E^ >P R) 



J AB 



c w 



5ba 



X B + e X A fi 5 ABX B i 



AB + & AB f L 



XBi 



-7 



-9f*xidjYZ» [D XY ]-ef x idJ Y A^ [Q?], 



where 



E 



E 



AB 

R(X~) 
AB 



E 



E 



L(x°) 

AB 

R(x°) 

AB 



-—((0 L ) A1 (Ol) B1 + l(0 L ) A2 (Ol) B2 ). 

-—((O r )ai(0* r ) b1 + \{O r ) A2 {O r ) B2 
cw V 2 

( t;(On)a3{0* n )b3 - t:(O n )a4(0* n ) B 4 

cw \2 2 

Mx )* 



—E 



AB 



D f x L Y = l( L (U f ) Xl (U* I ) Yl -Q f 3 



f~L S W 

cw 



(C.3) 
(C.4) 
(C.5) 

(C.6) 
(C.7) 

(C.8) 

(C.9) 
(CIO) 



Here, in definition of D xy index i is summed over generation indices 1, 2, 3. Q? L is a charge 
of the sfermion f^, and ig is its third component of isospin in the weak basis. The charges 
and third components of the isospin of the weak-basis fields for charginos and neutralinos 
are explicitly written in the definitions of constants E^'^ x ■* and E^'^ x ■* . 



D. Loop functions 

In this Appendix the loop functions appearing in the Z-boson amplitude and the box 
amplitudes are listed. 

Z-boson loop functions 

The Z-boson amplitude comprises two-loop functions from the triangle-diagram part of 
the amplitude, Fi(a,b,c) and F 2 (a,b,c), and two- loop functions from the self-energy part 
of the amplitude, fi(a,b) and f 2 (a,b) are 



Fi(a,b,c) = - 



1 



b-c 



a In a — b In b a In a — c In c 



^ , , v 3 11 
F 2 (a,b,c) = - - - 



46-c 



a — b a — c 

a 2 In a — b 2 In b a 2 In a — c 2 In c 



a — b 



a — c 



(D.l) 
(D.2) 
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1 In a -a 2 + & 2 + 2fe 2 (lna-ln6) 
fe ) = 2 " — + 4(^)2 ' (°- 3 ) 

, , 1 In 6 a 2 - 6 2 + 2a 2 (ln6- In a) 

AM) = 5 --5- + 4(^p ' < D - 4 > 

These loop functions are evaluated by neglecting the momenta of incoming and outgoing 
particles. The functions F\ and F 2 are symmetric regarding replacement of their arguments 
(a, b, c). In the limit of two equal argument the functions F\ and F2 have the following 
form, 

„ , , , s a — b — a In a + a In b /T ^ N 
F 1 (a,b,b) = ^— ^ , (D.5) 

„ , , 1 In 6 a 2 -6 2 + 2a 2 (ln6-lna) /n . 

F2(o , M) = ___ + _>_ 1. (D . 6) 

The arguments of logarithms appearing in the F\ and F2 can be divided by any constant, 
what can be used to redefine these functions as functions of two variables, for instance b/a 
and c/a. The unpleasant In 6 term in ji function can be replaced with ln(6/a) because of 
unitarity cancellations in the sum N^y^iAX* ^ ancj - therefore f\ and f 2 can be expressed 
in terms of one variable (b/a) only 



Box-loop functions 

Box amplitude contains two-loop functions, do and d 2 - 

xlnx ylny 
d {x,y,z,w) = - + 



(y - x)(z - x)(w - x) (x - y)(z - y)(w - y) 

zlnz w\nw 

(x — z) (y — z) (w — z) (x — w)(y — w )(z — w ) ' 

, .If x 2 In x y 2 In y 

d 2 {x,y,z,w) = - \ + 



4 I (y - x )i z - x )( w ~ x ) {x- y)( z - y)( w - y) 

^ z 2 lnz ^ w 2 In w 1 ^ 

(x — z)(y — z)(w — z) (x — w)(y — w)(z — w) J 

As mentioned before, they are also by evaluated neglecting the momenta of incoming and 
outgoing particles. 



E. Meson states and quark currents 

Meson states are assumed to contain valence quarks only. The quark-antiquark (q a <ll) 
content of the pseudoscalar-meson states is given in the table below. The quark-antiquark 
content of the vector- meson states is obtained replacing the fields K + , K , 7r + , 7r°, 7r~, 
K° K~, r/8, 771, 77 and 7/ by fields K* + , K*°, p + , p°, p~, ~K*® K*~ , fa, fa, (f> and u, 
and the angle 9p by the angle Qy. From the quark content of the meson fields one can 
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|M) 


quark content of \M) 


quark content of M(x) 


\K+) 
\K°) 
|vr+) 
|vr°) 

br-) 

1 / 

1%) 
\Vl) 


us c ~ bldt 

ds c 

ud c 

-^(uu c -dd c ) 

du c 

su c 

sd c 

-±(uu c + dd c - 2ss c ) 
-±(uu c + dd c + ss c ) 


su c ~ d s b u 

sd c 

du c 

-4(mz c - dd c ) 

ud c 

us c 

ds c 

-^(uu c + dd c - 2ss c ) 
-j^(uu c + dd c + ss c ) 


\v) 
W) 


cp\m) - sp\m) 
spivs) + c P \r)i) 


c P T] 8 (x) - sprn(x) 
sprjs(x) + cprn(x) 



Table 2: Quark content of the pseudoscalar-meson states and fields: The meson states listed in 
the Table 2 correspond to the tensor description of meson states The shorthand notation, 
cp = cos Op and sp = sin#p is used. 



find the meson content of e.g. axial vector (A) and vector (V) quark currents, (factor of 
proportionality, Lorentz indices and spinor structures are neglected), 



(5s> -(-7f-S)•' , + (-W + ^)•'' ,st^ ••' , + ^ , ™ 
^ ^ (t? " 70 *' + ( vf + 70 " tV ' s ^ + ^ + 

< I "-(-w-^)* ,+ (-w + ^) w, "^ + ^ <E - 2) 



Here sy = sin#y, and cy = cos 6y, and the numerical factors are normalizatons of mesons 
expressed in terms quark fields. The combinations of constants kg\ b are contained in ex- 
pressions for all quark currents (from the scalar-quark current to the tensor-quark current) 
and we introduce them to abbreviate the expressions for box form factors. 
The photon-penguin and the Z-boson-penguin amplitudes contain combinations of con- 
stants k?^ b given in Table 3. For instance, k% = Q u k^ u + Qdk p £, + Q s k§°. 
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k 


P° 






to 




kj 
k v z 


1 


%c 2W + 




7E SV 


Cy 


k 




V 










1 

V2 


Cp , Sp 




Sp Cp 

Vg 2V3 





Table 3. Combinations of constants appearing in photon-penguin and the Z-boson-penguin 
I t'P{V) amplitudes. 
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